Integration of natural data within a numerical model of 
ablative subduction: A possible interpretation for the 
Alpine dynamics of the Austroalpine crust. 

Manuel Roda 

Universita degli Studi di Milano, Dipartimento di Scienze della Terra "A. 

Desio", now at Universiteit Utrecht, Faculty of Geo sciences, Budapestlaan 4, 
3584 CD Utrecht, THE NETHERLANDS 

Corresponding author. E-mail: jediroda@gmail.com 

Maria lole Spalla 

Universita degli Studi di Milano, Dipartimento di Scienze della Terra "A. 
Desio", Sezione di Geologia, Via Mangiagalli 34, 20133 Milano - ITALY 
IDPA, CNR, Via Mangiagalli, 34, 1-20133 Milan, Italy 

Anna Maria Marotta 

Universita degli Studi di Milano, Dipartimento di Scienze della Terra "A. 
Desio", Sezione di Geofisica, Via L. Cicognara 7, 20129 Milano - ITALY 

IDPA, CNR, Via Mangiagalli, 34, 1-20133 Milan, Italy 

Short title: Model of ablative subduction in the Alps. 



1 ABSTRACT 

2 A numerical modelling approach is used to validate the physical and geo- 

3 logical reliability of the ablative subduction mechanism during Alpine con- 

4 vergence in order to interpret the tectonic and metamorphic evolution of an 

5 inner portion of the Alpine belt: the Austroalpine Domain. The model pre- 

6 dictions and the natural data for the Austroalpine of the Western Alps agree 

7 very well in terms of P-T peak conditions, relative chronology of peak and 

8 exhumation events, P-T-t paths, thermal gradients and the tectonic evolu- 

9 tion of the continental rocks. These findings suggest that a pre-coUisional 

10 evolution of this domain, with the burial of the continental rocks (induced 

11 by ablative subduction of the overriding Adria plate) and their exhumation 

12 (driven by an upwelling flow generated in a hydrated mantle wedge) could be 

13 a valid mechanism that reproduces the actual tectono-metamorphic config- 

14 uration of this part of the Alps. There is less agreement between the model 

15 predictions and the natural data for the Austroalpine of the Central-Eastern 

16 Alps. Based on the natural data available in the literature, a critical discus- 

17 sion of the other proposed mechanisms is presented, and additional geological 

18 factors that should be considered within the numerical model are suggested 

19 to improve the fitting to the numerical results; these factors include varia- 

20 tions in the continental and/or oceanic thickness, variation of the subduction 

21 rate and/or slab dip, the initial thermal state of the passive margin, the oc- 

22 currence of continental collision and an oblique convergence. 

23 Key words: ablative subduction; numerical modelling; hydrated mantle 

24 wedge; European Alps; Austroalpine Domain. 
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INTRODUCTION 

The exhumation of HP-U HP continental rocks is a topic that has intrigued 
many geoscientists since the late 1960 's. Based on depth-time paths analy- 
sis and analogue and numerical modelling, pre-coUisional and syn- to post- 



29 collisional exhumation mechanisms have been proposed (e.g. Duchene et al. 



1997). In the first case, exhumation occurs in an accretionary wedge (e.g. 



yamato et al. 2007), in regions of corner flow (e.g. Cloos, 1982; Shreve & 



Cloos, 1986 Cloos & Shreve, 1988) and in the serpentinised mantle chan- 



Gerya & Stockhert 2005 


Meda et al. 


2010 


Roda et al. 



2010). Crust-mantle delamination (e.g. Chemenda et al., 1995), continent or 



micro-continent collision (e.g. Gerya et al., 2008 Warren et al.. 2008), slab 



36 break-off (Ernst et al.. 1997), slab retreat (Ring & Layer, 2003), rollback 



iplab (Brun & Faccenna, 2008), intra-continental subduction (e.g. Schuster & 



Frank 1999 Thoni 2006; Janak et al. 2009 Stiiwe k Schuster 2010) and 



39 vertical extrusion leaded by rigid mantle indentation (RoUand et al., 2000 



Schwartz et al. . 2007; Schreiber et al.. 2010) are mechanisms that have been 



proposed to drive syn- to post-coUisional exhumation. 

Natural data and numerical predictions indicate that peculiar thermal 
and kinematic features characterise the uplift path of continental rocks ex- 
humed via different mechanisms. Rocks exhumed before the collision show 
cold prograde and retrograde paths, metamorphic peak imprints that are 
mainly recorded under HP and subordinate U HP conditions and contrast- 
ing metamorphic evolution recorded by adjacent portions of the same struc- 



tural nappe (Spalla et al., 1996 Gerya & Stockhert 2005 Engi et al.. 2001 



Meda et al. . 2010). Furthermore, pre-coUisional exhumation occurs within 
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low viscosity and/or low-density complexes, and the exhumation rate ranges 



51 from slow to moderate (Roda et al, 2010). 



Rocks exhumed during or after the onset of collision generally show 
warmer prograde and exhumation paths in which the retrograde trajectory 
can be isothermal or characterised by heating during decompression, depend- 
ing on whether the uplift occurs during early or late in the collision. In these 
cases amphibolite facies re-equilibration is mainly predicted with the local 
occurrence of partial melting. The metamorphic peaks record HP and U HP 



iponditions and the exhumation rates range from medium to high (Ernst & 



Liou 2008 Guillot et al. ^ 2009). 



Although these features have been widely described in HP-UHP oro- 
genic belts, some controversies persist regarding the individuation of the ef- 
fective exhumation mechanisms of subducted continental rocks. The debate 
regarding the mechanisms proposed for the exhumation of the Austroalpine 
continental crust in the European Alps is still open. The Austroalpine Do- 
main is mainly composed of continental metamorphic, igneous and sedimen- 
tary rocks derived from the distal part of the African passive continental mar- 
gin (the Adria plate), involved in the Alpine subduction since the Cretaceous 



(early Alpine, Spalla et al. (2010) and references therein). Although there is a 
general recognition that these rocks have been involved in Alpine subduction- 
collision, their contrasting tectono-thermal evolutions and structural setting 



ijiave generated different interpretations of Alpine geodynamics (e.g. Polino 



et al. . 1990; Spalla et al, 1996 Beltrando et al. , 2010), as discussed in the 



following section. The identification of the exhumation mechanisms affecting 
the orogenic belts requires the correct evaluation of the actual dimensions 



4 



of crustal and/or mantle portions that share a similar tectono-metamorphic 



evolution after their coupling (e.g. Spalla et al, 2010). This topic remains 



under debate within the Alpine (e.g. Spalla et al, 1996 Schwartz et al. 



2000 Rosenbaum & Lister, 2005) and extra- Alpine (e.g. Schneider et al. 



2004 Blanco- Quintero et al, 2011) scientific community. 



To contribute to this debate, a numerical model of ablative subduction 



Polino et al. 


1990 


Spalla et al, 


1996 


Meda et al, 


2010 


Roda et al. 



2010, 2011) with hydrated mantle wedge in an ocean-continent system is 
performed to compare, for the first time, the results with the natural data 
obtained for the entire Austroalpine Domain in terms of P-T peak estimates, 
relative peak chronology and exhumation time, P-T-t paths, thermal gra- 
dients and the geodynamic evolution of the continental rocks. In addition, 
we infer the volume of those parts that likely share the same tectonic and 
thermal history from peak to exhumation. Other proposed mechanisms are 
also discussed, based on the natural data collected and a brief comparison 
with present-day ocean-continent subduction in North America is proposed. 



91 THE AUSTROALPINE DOMAIN 

92 The uppermost tectonic domain of the Alpine belt is composed of Aus- 

93 troalpine continental rocks that tectonically overlie the Penninic Domain 

94 (Figs 1 and 2). The latter consists of mixed oceanic and continental slices 

95 and containing the remnants of the Tethyan domain. The western termi- 

96 nation of the Austroalpine units is located at the northern boundary of the 

97 Lanzo Ultramafic Massif and the eastern termination coincides with the east- 

98 ern border of the chain (Pannonian Basin). In the western portion of the 
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belt, the Austroalpine rocks also occur in small slivers mingled with conti- 
nental and oceanic units of the Penninic Domain. In the following sections, 
the tectono-metamorphic evolution of the Austroalpine Domain will be sum- 
marised by considering three different parts (Fig. 1): (1) the Austroalpine 
of the Western Alps outcropping to the west of the Lepontine Dome; (2) the 
Austroalpine of the Central Alps between the Lepontine Dome to the west 
and the Tauren Window to the east; (3) the Austroalpine of the Eastern Alps 
that outcrops to the east of the Tauern Window. 

The internal structure of the Austroalpine domain (Fig. 2), reveals a 
heterogeneous lithostratigraphy of the units, characterized by different meta- 
morphic evolutions, both in Western and Eastern sectors of the chain. The 
metamorphic evolution of the Austroalpine basement indicates that numer- 
ous portions, located along the entire length of the chain, have been deeply 



involved in subduction processes (e.g. Bousquet et al, 2004 Thoni, 2006) be- 
cause they display an eclogite- to ifP-amphibolite-facies metamorphism, as 



„well shown in the Map of the Metamorphic Structure of the Alps (Oberhansli 



[, 2004). The P/T ratios inferred for the eclogitic imprint generally increase 
from east to west, suggesting that the subduction setting of the Western 
Austroalpine is generally colder than the Eastern Austroalpine (Fig. 3a). 
However, highly-contrasting P/T ratios can be observed in the same struc- 
tural units in the entire Austroalpine Domain (i.e. the Sesia-Lanzo zone 
(SLZ), Fig. 3a). Metamorphic ages inferred for the eclogitic imprint of the 
Eastern Austroalpine are older than those obtained for the western part (e.g. 



Bousquet et al. 2004 Handy & Oberhansli, 2004), suggesting that the con 



123 tinental crust was involved in subduction in this portion of the Alpine belt 
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124 earlier than in the Western Alps (Fig. 3b). 



The Austroalpine domain of the Western Alps 

Alpine evolution in the Austroalpine domain of the Western Alps (AWA) is 
characterised by polyphase deformation under blueschist- to eclogite-facies 
conditions followed by retrogression under blueschist- to successive greenschist- 



facies conditions (e.g. Meda et al. 2010). For the Sesia-Lanzo zone (SLZ) 



mineral ages ranging between 60 and 80 Ma (Fig. 3 and Table 1) have been 



ijrelated to the Alpine eclogite facies peak (e.g. Lardeaux & Spalla, 1991 Inger 



IftaLl proi IDuchene et aTl \T997\ [Rubattol pMj [Rubatto efd] p^99| [Rebay 



& Spalla 2001 Rubatto et al, 2011). The P-peak values, which are widely 



134 dispersed, range between 1.0 and 2.5 GPa, and the maximum estimated value 



ipf T is approximately 650-680° (Fig. 4a). For the External Slices (Dal Piaz 



et al. , 2003 ) , lower P-peaks have been estimated with mineral ages ranging 



from 40 to 50 Ma (Fig. 3 and Table 1). 

Contrasting P-T and P-T-d-t evolutions have been inferred based on 



139 phase equilibria calculations and geothermometry (Castelli et al. (2007), and 



references therein) for different parts of the complex, which experienced a 



metamorphic imprint from eclogite to greenschists facies (e.g. Andreoli et al. 



19761 jPal Piaz et al.\ [Kienastj \T983\ [Tropper fc Essenej [2002| |Hellwig 



^003| IBabist et a/.j [2006| |Rebay fc Messigal [20071 [Roda fc ZucaHj [2008| [Gosso 



et al, 2010), Fig. 4). The P-T-t paths involve high dP/dT conditions and 



lawsonite can occur in assemblages developed during prograde or retrograde 



146 trajectories (Zucali & Spalla (2011), and references therein). The estimated 



ipxhumation rates lie in the range 1.4 to 3.5 mm/a (Zucali et al., 2002; Ru 



batto et al., 2011). Rubatto et al. (2011) suggest that multiple stages of 



re-equilibration occurred under eclogite facies conditions for the rocks of the 
Sesia-Lanzo zone, with two subduction-exhumation cycles occurring over 20 
Ma. Large-scale shear zones developed during the final stages of greenschist 



facies re-equilibration in central SLZ (Babist et al. 2006) and in the Dent 



iPlanche klippe (Pennacchioni & Guermani 1993 Hellwig, 2003 Roda & 



Zucali, 2008), and a later brittle- ductile faulting postdates the Oligocene 



andesitic dykes and igneous stocks of the Biella and Traversella emplace- 



156 ment in the SLZ (Zanoni et al, 2008 2010; Zanoni 2010). Both plutons are 



emplaced in the innermost part of the SLZ (the Eclogitic Micaschist Com- 
plex) close to the internal boundary of the Alpine metamorphic units (i.e., 
the Periadriatic tectonic line). 

Natural data related to the metamorphic peak conditions estimated for 
different portions of the Austroalpine of the Western Alps (reported in Figs 
la and 4a and Table 1) indicate the existence of a thermally-depressed en- 
vironment compatible with active oceanic subduction (the cold subduction 
of Cloos ( 1993[ )), with a low T/depth ratio also maintained during the ex- 
humation, which is testified by the blueschist re-equilibration following the 



eclogite facies imprint (e.g. Pognante, 1991). This peculiar metamorphic 



167 evolution suggests the interpretation of a burial-exhumation cycle that was 



ip,ccomplished during active oceanic subduction (e.g. Spalla et al. , 1996 Zucali 



et al, 2004 Meda et al. 2010; Roda, 2011). 



According to other authors, the occurrence of HP and UHP metamor- 
phism in the Western Alps can be interpreted as episodic and induced by 



fast changes from shortening to extension (Rosenbaum & Lister, 2005). In 



173 such a tectonic scenario, orogenic processes are governed by the subduction 
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174 of small ocean basins, slab rollback and the accretion of continental ribbons. 

175 In this tectonic configuration, each accretionary event is followed by revival 

176 of subduction and slab rollback in an adjacent oceanic basin, which, in turn, 

177 would switch the overriding plate into an extensional mode. The latter could 

178 then lead to the formation of metamorphic core complexes on the overriding 



179 plate and rapid exhumation of HP rocks (Rosenbaum & Lister (2005) and 

180 references therein). 

181 The Austroalpine domain of the Central Alps 

182 Peak P-T conditions for the eclogite facies event, in the Central Upper Aus- 

183 troalpine, have been derived from metabasics and their country rocks at 0.8~ 



, 1991 


Poh 


1991 


Spalla 


1993 


Solva 



et al. 2001; Habler et al, 2006). Zanchetta (2007) suggests the occurrence 



of UHP eclogites (2.6-2.9 GPa and 630-690°C) preserved as boudins within 
the garnet-amphibolites (Figs lb and 4b). Recent geochronological data 
for the age of Cretaceous metamorphism in the Upper Austroalpine nappes 



189 range between 76±5 and 85±5 Ma (Solva et al.. 2001 Habler et al.. 2006 



Zanchetta, 2007); an age of 143±2 Ma is inferred by Hoinkes et al. (1991) 
using Rb/Sr ratio on white mica (Fig. 3 and Table 2). A biotite-whole-rock 
Rb-Sr age of 70-80 Ma indicates that cooling below approximately 300°C 



193 occurred during the Upper Cretaceous (Habler et al. . 2006). 

194 During the Alpine time, a clockwise P-T path is inferred for the Err, 

195 Sella and Magna nappes: burial of the continental crust during the Ligure- 

196 Piemontese oceanic-crust subduction (between 100 and 80 Ma) and subse- 
quent multi-stage nappe stacking under high P /T greenschist facies condi- 



197 



198 tions are interpreted to have occurred between 80 and 30 Ma (Handy et al. 
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1996). Bachmann et al. (2009) indicate an age ranging from 54 to 47 Ma for 



the cessation of subduction induced by the Middle-Penninic (Brianzonese) 
colhsion. 



With the exception of the Languard-Campo ( Gazzola et al. 2000 ) and 



Err/Margna nappes (Handy et al.. 1996 Konzett & Hoinkes, 1996), the P~ 
T-t paths for this portion record a generally higher thermal state during 
the early stages of the exhumation with respect to the metamorphic peak 



206 (Hoinkes et al. , 1991 Habler et al. , 2001 Solva et al. , 2005 Zanchetta, 2007), 



reaching amphibolite facies conditions. 

The structural and metamorphic setting of this portion of the Austroalpine 
Domain has led to the interpretation of a continuous tectonic erosion of 
the upper plate in a subduction channel context regarding the Lower Aus- 



troalpine units (Bachmann et al., 2009). Conversely, for the Upper units, a 



syn-collisional extruding wedge is proposed by Solva et al. (2005). 



The Austroalpine domain of the Eastern Alps 

The P-T peak estimated for the Eastern Austroalpine exhibits a low P /T 
ratio, with values of T ranging from 550 to 780°C and P-peak values that 
vary widely from approximately 1.0 to 2.5 GPa (Figs lb and 4b). Excluding 
the Sieggraben eclogites and Reckner blueschists unit, there is a marked 
synchronism of ages for HP metamorphism along the Eastern Austroalpine 



of 80-90 Ma (Thoni (2006), Fig. 3b and Table 2). The slight difference with 



respect to the peak-age of the Western Alps led Thoni (2006) to infer that the 



HP event in the Eastern Alps resulted from a distinct event of continental 



ipubduction, probably related to the Meliata-Hallstatt ocean closure (Kurz & 



Fritz 2003 Thoni 2006; Handy et al. 2010). However, Handy et al. (2010) 
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note that the Eo-alpine orogenic wedge contains no suture in the classical 
sense of an ophiolite belt sandwiched between metamorphosed upper and 
lower-plate units. The only relics of the Meliata-Hallstatt ocean consist of 



ipelagic sediments incorporated in Jurassic-age accretionary wedges (Mandl 



& Ondrejicka, 1993). 



The lack of sedimentary material sourced from high-grade metamorphic 
rocks in the Late Cretaceous (Campanian - Maastrichtian) central Alpine 
Gosau sediments may indicate a decrease of exhumation rates after approx- 
imately 85 Ma (Wiesinger et al. , 2005; Thoni[ 2006) which has been inter- 
preted as the result of a short-lived (5-7 Ma) high-P event, followed by rapid 



234 exhumation (Janak et al., 2009). This view is partly contradicted by that 



235 of Tenczer & Stuwe (2003), who show that northern Koralpe and the sur- 



rounding Austroalpine crystalline basement cooled to below 200° C at 50 Ma 
(fission track), suggesting a lower exhumation rate than those invoked above. 
The P-T-t paths for this portion are characterised by a heating during 



239 the exhumation or by an isothermal uplift. Alternatively, Gaidies et al 



240 (2008a) report an anti-clockwise Alpine P-T path in the Rappold Complex, 



demonstrating a thermal state during exhumation that was colder than that 
of the burial event. 

The exhumation of Eastern Austroalpine units could be related to the 
occurrence of a back-arc extension during the Ligure-Piemontese subduction 



245 successive to the Meliata closure (Kurz & Fritz, 2003) or driven by an ex- 



246 truding wedge and erosion during a transpressive regime within an oblique 



247 subduction event (Thoni, 2006; Handy et al. . 2010). 
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248 NUMERICAL SIMULATION 



249 Model setup 



250 Submar code (Marotta et al, 2006) is used to simulate the convergence of 

251 an ocean-continent system using a grid 1400 km wide by 708 km deep con- 

252 taining 6165 nodes and 3010 triangular quadratic elements as described in 



Roda 


(2011 


) and 


Roda et al. 


(2011 



254 laws are summarised in Table 3. Crust and mantle are compositionally dif- 



255 ferentiated via the Lagrangian particle technique ( Christensen , 1992 ) with a 

256 density of 1 marker (particle) per 0.25 km^. To simulate plate convergence, 

257 the velocities are fixed at the base of the oceanic crust (a depth of 10 km) 

258 with a value of 3 cm/a. To force subduction, we fixed the same velocity at 

259 depths up to 100 km at the nodes of the numerical grid distributed along a 

260 45°dip plane. We assumed a continental lithosphere thickness of 80 km (with 

261 30 km of continental crust) to represent an originally-thinned passive margin 



that characterised the former margin of Adria (e.g. Dal Piaz, 2001). Using a 
value of approximately 40 Ma of active spreading for the Ligure-Piemontese 



ocean (Handy et ai, 2010), we assumed 80 km thick oceanic plate following 



265 the relation ( Turcotte & Schubert , 2002 ) : 



y = 2.32(H) 



1/2 



(1) 



266 where y is the thickness of the oceanic lithosphere, k is the thermal dif- 

267 fusivity (10~^ iv? / s) and t is the oceanic plate age. 

268 The model runs for 65 Ma of active subduction at time increments of 

269 100 ka; the model simulation ends before continental collision. The assumed 
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model constraints cause a strong coupling (i.e. high friction) between the 



upper and lower plates (see Roda et al. (2010) for a discussion of this and 



272 additional references); therefore, our models are representative of the sole 

273 erosive margin type, which represent, at present, approximately the 60% of 



the global active margins (Clift & Vannucchi, 2004). 



275 Model results 

276 Figure 5 illustrates the predicted temperature and velocity fields, and the 

277 distribution of markers at four time steps. In the early stages of the sub- 

278 duction the ablation of the overriding plate by the oceanic plate dominates 

279 the subduction dynamics. The high coupling between the two plates inhibits 

280 the formation of a large and deep trench, and only a small amount of sed- 

281 iment is generated (Fig. 5a). The kinematics are dominated by the burial 

282 flow, although some small-scale upwelling flows could develop in the mantle 

283 wedge (Fig. 5a, inset). However, the intensity of the upwellings are not able 

284 to exhume a large amount of subducted material. With the progression of 

285 mantle wedge serpentinisation (which is due to the continuous subduction 

286 and dehydration of the oceanic plate, up to 20 Ma from the beginning of 

287 subduction), a more intense upwelling flow develops, especially in the upper 

288 part of the mantle wedge (Fig. 5b). Consequently, some crustal and sedi- 

289 mentary material is carried up to superflcial levels. The sub-horizontal flow 

290 generated along the overriding plate maintains the high degree of coupling 

291 between the two plates and induces a flexural bending of the continental 

292 plate close to the trench with consequent uprising of the lithospheric mantle. 

293 The suction generated by the high degree of coupling decreases the slab dip 

294 at depths of 200-350 km (Fig. 5b). During the intermediate stages of sub- 
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295 duction the upwelling flow in the serpentinised mantle wedge becomes the 

296 dominant mechanism, the coupling between the two plates decreases and a 

297 large trench wedge develops (Fig. 5c). Large amounts of subducted mate- 

298 rials are exhumed, and the upper and lower continental slices and sediment 

299 particles are involved in the recycling of the mantle wedge. However, only 

300 part of the recycled material reaches the surface; the rest remains in the deep 

301 portion of the mantle wedge or is more deeply buried (Fig. 5c, light purple, 

302 deep purple and yellow markers). Due to the decrease in plate coupling, suc- 

303 tion decreases and the slab dip increases in comparison to the previous time 

304 increments. The upwelling flow increases the rising of the lithospheric mantle 

305 to subsurface levels, and thermal and mechanical erosion of the continental 

306 plate occurs. 

307 Kinematically, the mantle wedge is characterised by three vortical flows 

308 located at depths from 30 to 100 km (Fig. 5c, inset). The late stages are 

309 characterised by a stable and highly efficient upwelling flow arising from 

310 the widespread serpentinisation of the suprasubductive mantle wedge. The 

311 consequence is a large orogenic wedge generated by the accretion of exhumed 

312 upper and lower continental crustal slices, recycled trench sediments and 

313 particles of upper oceanic crust. Even markers of continental lithospheric 

314 mantle are buried and exhumed in the mantle wedge (Fig. 5d, see markers 

315 legend). The combined effect of ablation and upweUing flow induces the 

316 thermal and mechanical erosion of the overriding plate along a zone extending 

317 from 100 to approximately 150 km from the trench, with the consequent 

318 denudation of the lower crust (Fig. 5d). The streamline pattern shows a 

319 single and intense vortical flow located in the mantle wedge from a depth of 
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320 30 to approximately 100 km (Fig. 5d, inset). 

321 The P-T-t paths followed by the upper and lower continental markers 

322 indicate burial of the material during the early stages of subduction and sub- 

323 sequent exhumation to subsurface levels and the upper part of the mantle 

324 wedge (Fig. 6). The mechanism of burial depends on subduction age; the 

325 T/P ratio decreases during the intermediate stages depending on the coohng 

326 generated by the progressive subduction of the oceanic plate. Furthermore, 

327 because the trajectory of every marker could be characterised by more than 

328 one burial-exhumation loop as suggested by the P-T-t paths (Fig. 6), the 

329 markers could record different values of Prmax a-^d Tj^^x during their evolu- 

330 tion within the mantle wedge and different exhumation rates. 

331 Several analyses can be performed taking the exhumed markers into ac- 

332 count; in this model, all particles that reach a minimum burial depth of 80 km 

333 and are subsequently exhumed to a depth of at least 30 km are considered as 

334 exhumed markers. Analysis of the thermo-barometric conditions recorded by 

335 the exhumed markers confirms the high variability of the thermo-barometric 

336 conditions at maximum pressure (PmaxT), with Pmax values from 1.2 GPa to 

337 up to 3.7 GPa, and of T, ranging from 250 to 600°C (Fig. 7a). The distri- 

338 bution of the predicted thermo-barometric conditions recorded at maximum 

339 temperature by the markers {PTmax, Fig- 7b) also demonstrates high vari- 

340 ability of the pressure associated with Tj^ax (from 1.0 to 3.5 GPa), but less 

341 variation in temperature, which is limited to the range 450-650°C. At the 

342 end of the simulation, the upper 50 km of the accretion wedge is composed 

343 of 56% mantle and hydrated mantle, 43.5% continental material and 0.5% 

344 oceanic material. Markers belonging to the lower continental crust are about 
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345 17.5% of the total amount of the exhumed material whilst 10% represents 

346 the amount of upper continental crust. 

347 The thermal-gradients characterising the subduction zone, which are pre- 

348 dieted at various reference depths from the minimum values of the depressed 

349 geotherms, (Fig. 7c), exhibit the strongest decrease during the early stages of 

350 their evolution, followed by a modest decrease during the intermediate stages 

351 until a constant value is reached during the late stages. If we consider the 

352 sole mantle wedge area (up to 300 km depth), the thermal gradients during 

353 the late stages of subduction range from 2 to 4 °C/km. Thermometric deter- 



354 mination of UHP orogenic peridotites (Ernst & Liou, 2008) corroborate such 

355 predicted cold thermal regimes. Furthermore, the mean thermal gradients 

356 obtained from the burial path of the markers in the numerical simulation 

357 overestimate those extracted from the isotherms, at least during the first 40 

358 Ma of subduction (Fig. 7c). Analysis of the P/T ratio, which is obtained 

359 from the mean PmaxT conditions recorded by the exhumed markers at five 

360 different time steps (Figs 7d and e), confirms the progressive cooling of the 

361 mantle wedge over time during ongoing subduction; the P/T values, which 

362 range from approximately 2.5 to 3.6-3.7 (10^^ GPa/°C), are equivalent to 

363 T/depth values ranging from 13 to 9 °C/km. It is important to emphasise 

364 that these P/T ratios do not coincide with the thermal gradients, which 

365 are exclusively represented by the configuration of the geotherms at different 

366 time steps (Fig. 7c). 

367 We calculated the dimensions of the crustal markers groups that expe- 

368 rienced common P-T-t histories, from the peak to the exhumation (60 Ma 

369 from the beginning of the subduction), disregarding their origin (upper or 
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370 lower crust). Taking into account that each marker represents a homoge- 

371 neous area of 0.25 km^, we found a maximum group dimension of 53 km^ 

372 with a mean value of 11 km^. These estimates represent the dimensions of a 

373 single metamorphic unit that means a portion of crustal and/or mantle ma- 

374 terial that shared a common tectono-thermal history within the exhumation 

375 process. 

376 COMPARISON WITH NATURAL DATA 

377 In this section a comparison between model predictions and natural data from 

378 the Austroalpine continental crust is presented and the results are critically 

379 discussed. Natural data are grouped according to the locality and tectonic 

380 units, P-T estimates and peak ages (Figs 1 and 4, Tables 1 and 2). We 

381 consider all the particles that have been buried below and recycled above 

382 a depth of 30 km (1.0 GPa), which is the starting crustal thickness, to be 

383 exhumed continental crust. For this reason, no statistical data are presented 

384 below a pressure of 1.0 GPa. 

385 The Austroalpine of the Western Alps 

386 In Fig. 8, the thermo-barometric peak conditions of the natural data for the 

387 Austroalpine of the Western Alps are plotted against the numerical results 

388 in terms of of Pressure- (Fig. 8a, b) and Temperature- (Fig. 8c, d) peak 

389 conditions. Although there is in general good agreement between the P-peak 

390 recorded by the rocks and that of the model, the better agreement in terms 

391 of temperature is obtained by comparison with the numerical P-T peak 

392 conditions referred to the maximum temperature {PTmax)- The predicted 

393 values of T at Pmax are generally lower than those derived from the natural 
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394 assemblages, particularly for T-values related to higher maximum pressures 

395 (from 1.5 to 2.5 GPa), whereas the T-values estimated for lower maximum 

396 pressures (<1.5 GPa) fit well. The comparison of natural data with model 

397 results in regard to the PT^ax are instead ever in very good agreement. 

398 This agreement could indicate that the P-peak recorded by the Austroalpine 

399 rocks does not correspond to the maximum pressure reached during the rock 

400 burial, but could instead be related to the pressure conditions recorded at 

401 the maximum temperature, where the kinetics of the metamorphic reactions 

402 increase and the phase changes can be more pervasive. Conversely, this could 

403 also be due to the lack of sufficiently refined petrological methods able to 

404 identify actual pressures instead of minimal pressure conditions and reduce 

405 the uncertainty in estimates of T (Fig. 4). 

406 If the subduction simulation has a duration that is equivalent to a time- 

407 span of 65 Ma, i.e., starting at 95 Ma and ending at 30 Ma, it is possible 

408 to compare the simulation results with the timing of P-T peak indicated in 

409 the literature. These ages range from 80 to 60 Ma for the SLZ and from 

410 50 to 40 Ma for the External Slices (see Table 1 and Figs 2, 8e and f). 

411 Although older peak ages were proposed in the literature, the improvement 

412 of the dating techniques and the problem related to the excess and/or loss 



413 of Argon in the Alpine belt (e.g. Ruffet et al, 1995) suggest that the more 

414 recent estimates available and not related to Ar methods are more reliable 

415 (see Table 1). Based on geological evidence, part of the subducted material 

416 was already exhumed to shallow crustal levels by 30 Ma, and this result is 

417 compatible with the development of greenschist facies assemblages at the 

418 expense of eclogite facies ones when the subduction ceased at the beginning 
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419 of continental collision (e.g. Piatt 1986 Polino et al., 1990 Zanoni et al. 



2010). Following these natural contraints, the end of the model simulation 



421 (after 65 Ma, i.e., 30 Ma) is considered to broadly coincide with the age 

422 at which exhumation occurs. Based on these assumptions, there is very 

423 good agreement temporally between the simulated P-T conditions and the 

424 record of the peak conditions in the natural rocks for both the SLZ and 

425 External Slices (Figs 8e and f). At the end of the simulation, i.e., at 30 Ma, 

426 a large number of the same continental ablated and buried markers have been 

427 exhumed to subsurface levels (Figs Be and f, yellow markers). Some of the 

428 markers, however, reside at significant depth at the end of the simulation, and 

429 could potentially remain at depth unless exhumed by a subsequent orogenic 

430 event. 

431 In Fig. 9, a comparison between selected natural and simulated P-T-t 

432 paths is proposed that also takes pre- Alpine P-T paths into account, where 

433 available in the literature. This is useful for constraining the extraction of 

434 the marker trajectories with respect to their starting structural level, at the 

435 beginning of the subduction and in the upper or lower crust, based on the 

436 results suggested by the pre-Alpine P-T-t natural evolution. Furthermore, 

437 we extract the simulated P-T-t paths by considering the time interval cor- 

438 responding to the specific Alpine peak age proposed in the literature or, in 

439 the absence of radiometric data regarding natural peak assemblages, the P- 

440 T-t path recorded during a 80 to 60 Ma time span. Although many P-T-t 

441 paths are reported in the literature, we compared the numerical results with 

442 P-T-t paths produced since 2000 for which modern and accurate structural 

443 and petrological analyses were available. In Fig. 9 the grey area represents 
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444 the area comprising all the numerical P-T-t trajectories, and the purple and 

445 blue lines represent some examples of various paths. Orange and cyan boxes 

446 and black arrows refer to the natural P-T conditions (pre- Alpine and Alpine) 

447 and natural P-T-t paths, respectively, as given in the literature and noted in 

448 the figure legend. Natural and simulated P-T-t paths are generally in good 

449 agreement in terms of both prograde and retrograde trajectories with the ex- 

450 ception of the External Slice paths Arolla and Mt. Morion (Fig. 9). In these 

451 cases, as discussed in the literature, peak conditions are attained close to 



452 the onset of continental collision at approximately 30 Ma, (e.g. Zanoni et al. 



2010 Handy et al. . 2010) and consequently the thermal regime, character- 



454 ising the exhumation of these rocks, can increase due to high shear heating 



455 and strong radioactive decay (e.g. Gerya et al., 2008 Warren et al, 2008). 



456 It must also be emphasized that our numerical simulation is stopped before 

457 the collision event, and therefore, exhumation is only partially accomplished 

458 under a sensible low-temperature regime. 

459 The Austroalpine of the Central Alps 

460 The comparison between the natural and the simulated thermo-barometric 

461 conditions of the Central Alps shows no agreement with regard to PmaxT 

462 (Figs 10a and b) values because of the higher temperatures retrieved from 

463 the natural data. However, the natural data are consistent with the higher 

464 simulated T^ax values in the PT^ax diagram (Figs 10c and d). The samples 

465 that record lower peak- pressure (lower than 1.0 GPa) are not compared to 

466 the simulations. 

467 For the comparison between the natural age data and the timing of the P- 

468 T peak from the model, we extract the simulated P-T values for the first 20 

20 



469 Ma of subduction (i.e., from 100 to 80 Ma, as suggested in the literature. Fig. 



iPb and Table 2), and the exhumation age corresponds to 50 to 35 Ma (Handy 



et al. , 1996; Bachmann et al. , 2009). Fig. lOe allows a qualitative estimate of 



fitting percentage between simulated P-T peak values (red dots) and natural 
estimates (shaded area) of at least 50%. Furthermore, the predicted P-T 
conditions at the end of simulation (orange dots) are coherent with natural 
P-T data set and associate age. 

Only one natural P-T-t path agrees well with the predictions of the 
model (Fig. 9: "Mortirolo" (ref. 1 Central Alps, Fig. lb)). Other P 
T paths from the literature display a warmer thermal state for both the 
peak and exhumation conditions and are therefore not consistent with the 
simulated P-T-t paths. 

The Austroalpine of the Eastern Alps 

A comparison of the natural and the simulated thermo-barometric conditions 
of the Eastern Alps shows no agreement with regard to the PmaxT values 
and a poor agreement with PTmax values due to the higher thermal state 
suggested by the natural data (Fig. 11). For the Eastern Alps, we considered 
130 Ma to be the age at which subduction began, based on the peak age of 
Sieggraben (Table 2). In this region, in contrast to the the Western Alps, 
the P-peak ages lie in the range from 100 Ma to 90 Ma, and two exhumation 



iptages are considered: 80 Ma (Thoni et al.. 2008) and 50 Ma (Tenczer & 



Stuwe, 2003). 



A low correspondence is obtained between the timing of the P-peak for 
the Eastern Austroalpine (90-100 Ma) and the model predictions (Fig. lie). 
Regarding the peak age interval, good agreement is obtained only for lowest 
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494 P values. For the higher pressures, a colder thermal state is predicted by 

495 the simulation (Fig. lie). A significant better agreement is found if the 

496 subduction beginning is assumed at 100 Ma (excluding Sieggraben), and the 

497 predicted P-T conditions refer to the early stages of subduction (1-10 Ma 

498 = 100-90 Ma, Fig. llf). Indeed, for this time span, a higher thermal state 

499 is predicted than for the intermediate stages of the subduction. In any case, 

500 a considerable amount of buried material has already been exhumed at both 

501 exhumation time-steps considered, although a large part sinks to greater 

502 depths and probably will remain buried. 

503 As with the Central Austroalpine, only one natural P-T-t path agrees 

504 well with the predictions of the model (Fig. 9: "Sopron" (ref. 19.2 Eastern 

505 Alps, Fig. Ic)). In this case, only the peak-pressure is well reproduced, 

506 whereas the thermal state of the exhumation path is underestimated by the 

507 simulated paths. Other natural P-T paths display a warmer thermal state 

508 for both the peak and exhumation conditions and are not reproduced by the 

509 simulations. 

510 DISCUSSION 

511 The model predictions agree well with the natural P-T estimates and the 

512 P-T-t paths derived from for the Austroalpine crust of the Western Alps. 

513 All natural P-T data are reproduced by the numerical P-T predictions, es- 

514 pecially the P-T peak conditions obtained by PTmax analysis. Even the time 

515 interval between the peak and exhumation conditions are closely consistent 

516 with the ages inferred from the rocks, with different P-peak ages, for the SLZ 

517 and the External Slices. In this context, the exhumation rate suggested by 
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518 the numerical model ranges from 0.1 to 1 cm/a, with mean values of 0.3-0.4 



cm/a. These exhumation rates agree with those predicted by Roda et al. 



520 (2010) and suggested by Zucali et al. (2002) and Rubatto et al. (2011) for a 



part of the SLZ. Good agreement is found between the natural and simulated 
P-T-t paths, both with and without consideration of their pre- Alpine evolu- 
tion, with few exceptions. These cases are characterised by a peak-age very 
close to the collisional event, and their retrograde paths could be affected by 
the increase of temperature induced by the continental collision. 

The same heterogeneity that characterises both the natural and predicted 
P-T peak conditions and the P-T-t evolutions in the same structural domain 
suggests the juxtaposition of rocks with different metamorphic evolutions. 
The rocks have been coupled and decoupled at different stages during their 



530 trajectories in the mantle wedge (e.g. Spalla et al., 1996; Gerya & Stockhert 



2005 


Bousquet 


2008 


Roda et al. 2010 


Spalla et al. 


2010) 



is evident that the tectono-metamorphic unit definition requires an evaluation 
of the actual size of crust and mantle-derived slices that coherently shared, 
after amalgamation, the same tectonic and thermal history. In light of this, 
the structural and metamorphic evolutions of basement rocks trace their 
transit throughout different levels of the lithosphere and sub-lithospheric 



ipantle (e.g. Schwartz et al, 2000 2007 Spalla et al, 2010 Blanco-Quintero 
et al.\ [2011 ), rather than the purely litho logic associations or the common 



dominant metamorphic imprint. The tectono-metamorphic unit dimensions 
obtained in this work are at least one order of magnitude lower than those 



proposed for the Alps (e.g. Oberhansli , 2004): this means that the refinement 



542 of a method for the identification of contours of thermally-characterized and 
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543 structurally-distinct units became crucial (Spalla et al. , 2005, 2010) to obtain 

544 a good definition of the Alpine tectono-metamorphic units. In the case of the 

545 Western Austroalpine, the model results and their fitting with the natural 

546 data indicate that a nappe (Sesia-Lanzo Zone) can result from the accretion 

547 of different tectonic units. 

548 The thermal gradient of the simulated subduction zone lies in a range 



549 from 2 to 4 °C/km (Fig. 7c). [Groppo et al\ ( |2009D , based on the P-T as- 

550 semblages obtained for some samples belonging to the "Lago di Cignana" 

551 unit, suggest a gradual decrease in the thermal gradient (from 9-10 to 5-6 

552 °C/km) over time. This variation, which is simply related to the time, has 

553 been interpreted by these authors as evidence for an increase in the subduc- 

554 tion rate of the Ligure-Piemontese oceanic slab during the Eocene. However, 

555 a deeper analysis of P-T-t paths reveals that the attribution of thermal gra- 

556 dient significance to a single and transient P-T estimate cannot be directly 

557 related to the variation of the geotherms. The inference of variation in ther- 

558 mal gradients requires a minimum of two P-T estimates and the knowledge 

559 of the configuration of isotherms in an evolving dynamic system such as an 



560 active subduction zone. Thus, the prograde path inferred by Groppo et al. 



561 (2009) suggests a thermal gradient of approximately 3 °C/km, which agrees 

562 well with that predicted by the numerical model, rather than the transition 

563 between two different thermal gradients proposed by the authors. Further- 

564 more, the numerical model suggests that the thermal gradients inferred by 

565 the burial path of the markers are quite different than those extracted from 

566 the isotherms (Fig. 7c). This finding indicates that the determination of 

567 thermal gradients that are derived directly from the P-T path could be af- 
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568 fected by uncertainties arising from the nonlinear trajectories of the markers 

569 in the mantle wedge and to variability in the thermal regime. 

570 The model results suggest the occurrence of a very dynamic mantle wedge 

571 that produces two or three subduction-exhumation loops recorded by the 

572 markers within 65 Ma of active subduction (Fig. 6). This aspect of the mod- 

573 elling is consistent with the recent discovery of two subduction-exhumation 

574 cycles recorded in the same rocks in the Western Austroalpine of the Alps (the 



575 Sesia-Lanzo Zone (Rubatto et ai, 2011). Similar subduction-exhumation cy- 

576 cles have been invoked to justify the complex P-T trajectory inferred for the 



577 rocks of the Caribbean subduction channel ( Blanco- Quintero et ai, 2011) 



and the Franciscan complexes (Wakabayashi & Dumitru, 2007 Wakabayashi 



2011). 



580 All of the previous points suggest a pre-collisional evolution of the Aus- 

581 troalpine of the Western Alps, involving the burial of slices ablated by the 

582 margin of the overriding continental plate (the Adria plate) during subduc- 

583 tion and their successive exhumation (driven by the upwelling flow generated 

584 in a hydrated mantle wedge). 

585 A much more modest agreement between the natural data and model 

586 predictions occurs for the evolution of Austroalpine rocks from the Central 

587 and Eastern Alps. In these regions, the simulated thermal state is generally 

588 underestimated with respect to that derived from P-T estimates of the rocks. 

589 Although the timing of simulated P-T peak and exhumation are comparable 

590 with those suggested in the literature, only a small part of the predicted P-T 

591 conditions is consistent with the natural conditions. 

592 These results suggest that entirely pre-collisional ablative mechanisms are 
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unable to explain the high temperature recorded by the rocks belonging to 
the Austroalpine of the Eastern Alps. A model that includes continental col- 
lision could better explain the high thermal state reached by the rocks before 



2008 


Gerya et al., 


2008 


Warren 



et al. , 2008 ), although the high levels of mixing of units predicted by the cited 



598 models do not appear to reproduce the low level of heterogeneity observed in 



599 the Eastern Alps, according to Janak et al. (2009). An intra-continental sub 



duction zone is the suggested alternative interpretation (Janak et al.. 2009). 
However, a mechanism able to engage an intra-continental subduction is dif- 
ficult to identify. The thickening of the continental crust of the upper plate 
is probably a more common mechanism (which is effective during active sub- 
duction) that could be invoked to explain higher values of Tmax- Indeed, the 
influence of a thick upper plate on the thermal state of the wedge has been 
described in the case of low slab dip and strong plate coupling, where a higher 
thermal regime characterises the mantle wedge during active subduction be- 



fore continental collision (Roda et al. . 2011). This flnding suggests that the 



higher T values recorded in the Eastern Austroalpine could be related to lat- 
eral variations of upper plate thickness along the Alpine convergent margin. 
As noted previously, the thermal state depends on the subduction age (early 
vs. intermediate subduction. Figs 7d and 11). 

The poor lithologic heterogeneity assessed in the Eastern Alps is weak- 
ened by the heterogeneities in the P-T peak conditions and P-T-t evolutions 
of natural rocks both within single and between different units (Fig. 1 and 
Table 2). The lithostratigraphic homogeneity that characterises the East- 
ern Austroalpine with respect to the melange-like assemblage of the Western 
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Austroalpine does not necessarily imply homogeneous and coherent P-T-t- 
d paths. The lithostratigraphy of the Western Austroalpine is comparable 
to that described in the Central-Eastern Austroalpine, but a high level of 
heterogeneity is predicted by the different P-T-t-d evolution recorded by 
the slices constituting the same structural domain (or "Alpine nappe") (e.g. 



Schwartz et ai, 2000 Bousquet, 2008 Spalla et ai, 2010). In this con 



text, better-refined prograde P-T-t paths and, in general, a wide spread 
in derived P-T-t or P-T-t-d paths are needed for the Austroalpine of the 
Central-Eastern Alps to allow the existence of possible tectono-metamorphic 
homogeneities to be discerned. 



Roda et al. (2010, 2011) suggest that a thinned oceanic plate (less than 



80 km, as reasonable for a young ocean) and/or a variation of the subduction 
rate and/or the slab dip could affected the thermal state of the mantle wedge. 
Taking into account an oblique Alpine subduction and the possible V-shape of 



632 the Ligure-Piemontese ocean (e.g. Schmid et al, 2004), it should be possible 



to obtain a decrease of the thermal state of the subduction zone from east 
to west that is related to an increase in the oceanic plate thickness. 

However, some natural P-T conditions reflect a very high thermal state, 
partially consistent with the numerical predictions obtained for the early 
stages of the subduction. This observation suggests that the initial LP 
thermal state of the rocks could influence the thermal evolution of buried- 



639 exhumed rocks during the subduction. Marotta et al. (2009) (and references 



therein) indicate that the Austroalpine domain is widely affected by HT-LP 
Permian- Triassic metamorphism, which developed during the early stages 
of the Tethyan rifting. The subduction of the Eastern Austroalpine crust 



27 



643 earlier than that of the Western Austroalpine crust could imply a thermally 

644 perturbed continental margin, a younger (and therefore hotter) oceanic plate 

645 and a higher T / P ratio in the early stages of the subduction. 

646 Moreover, horizontal movement linked to transpressive subduction could 

647 affect the P-T conditions recorded during the exhumation paths, as noted 

648 by Thompson et al. ( 1997a||b ). In this case, the exhumation is characterized 

649 by a transition from blueschist to granulite fades conditions consequent to 

650 the increase of the transpressive component. 

651 For the comparison between natural data from the Austroalpine domain 

652 of the Alps and the model predictions, we briefly compare the simulation 

653 results with the geodynamic scenarios proposed for still active subduction 

654 settings such as Cascades and Franciscan. In these contexts, the continen- 

655 tal collision does not affect the subduction dynamics as in the Alps. The 

656 thermal structure of the Cascades is not so different with respect to the sim- 



657 ulated scenario: indeed, if we refer to fig.2 of Hyndman et al. (2005), notable 

658 similarities between the two models can be found. In the Cascades a cold 

659 fore-arc and thermal erosion of the continental lithosphere beneath the back- 

660 arc region are present. The same thermal structure is also predicted by our 

661 model, with a cold fore-arc region almost 100 km wide and the rising of the 

662 1400 K isotherm up to 50 km depth beneath the back-arc region, accounting 

663 for a high thermal regime. As a consequence, we do not expect differences 

664 concerning the exhumation mechanisms but a warmer exhumation-related 

665 P-T evolutions would be consequent to the shallower position of the 1600K 

666 isotherm (asthenosphere-lithosphere boundary). Another major difference 

667 between the two models is the larger extension of the back-arc region in the 
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668 Cascades. This could be related to the longer subduction activity of the 



Cascades (180-160 Ma) with respect to 65 Ma in our model (Currie et al. 



2008), that implies a continuous thermal erosion of the continental upper 



plate. Furthermore, the higher subduction dip of the Cascades (medium and 



672 deep dips: Lallemand et al., 2005 Cruciani et al.. 2005 Currie et al.\ 2008) 



673 with respect to the presented model, would induce a larger back-arc hot re- 



igion due to a very efficient convective cell beneath the upper plate (Currie 



et al. 


2008 


Roda et al. 


2010 



676 the Cascades scenario (Hyndeman, 1995; Flueh et al. . 1998) also implies a 



different thermal state with respect to that obtained by our model. Finally, 
the oblique subduction that characterizes large part of the Cascades could 



^generate a perturbed thermal regime within the subduction zone (Thompson 
et a/.[ |1997a||b| [Rondenay et q/.[|2008| ). 



Although an early intra-oceanic subduction scenario is invoked for the 
Franciscan compared with the ocean-continent setting adopted in our simu- 
lations, a lot of similarities between natural and simulated data are recogniz- 
able: (i) a wide range of peak metamorphic ages recorded by the exhumed 



rocks (Wakabayashi & Dumitru, 2007); (ii) different P-T evolutions resulting 



from different tectonic trajectories during subduction and exhumation, with 
the consequent accretion of different tectono-metamorphic units at different 



times in the same subduction zone (Wakabayashi & Dumitru, 2007); (iii) 
the record of multiple (2 or 3) burial-exhumation cycles of blue-schist rocks 



690 within 65 Ma of active subduction (Wakabayashi, 2011 ), although these rocks 



691 have been related in the past to an accretionary wedge / subduction chan- 



ipel mechanism rather than to a hydrated mantle wedge dynamics (Cloos & 
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Shreve, 1988); (iv) the presence of HP-HT metamorphic rocks, which show 



694 anti-clockwise P-T-t paths hkely related to the earliest steps of the subduc- 



695 tion zone (Smart & Wakabayashi , 2009) that are coherent with the predicted 

696 P-T-t of Fig. 6; (v) prograde eclogitization of a relatively cold subducting 

697 slab, and subsequent exhumation and blueschist-facies recrystallization of 

698 HP rocks characterized by exhumation rates of several mm/a, as suggested 



699 by the numerical results (Roda et ai, 2010; Roda, 2011). 



700 CONCLUSIONS 

701 The general very good agreement between the predictions of the model and 

702 the natural data obtained for the Austroalpine of the Western Alps suggests 

703 that the ablation of the upper plate lithosphere coupled with the dynamics 

704 of the hydrated mantle wedge is a valid pre-collisional mechanism that can 

705 reproduce the natural burial and exhumation evolution of this portion of the 

706 Austroalpine domain during active subduction. Although this model appears 

707 to be less compatible with the natural data of the Central-Eastern Alps, 

708 a variation of the oceanic and/or upper plate thickness along the Ligure- 

709 Piemontese ocean, a variation of the subduction rate and/or the slab dip, 

710 the initial thermal state of the passive margin, the occurrence of continental 

711 collision or an oblique subduction could justify a variation in the thermal 

712 state from east to west along the Alpine convergent margin. Furthermore, 

713 the simulated ocean-continent dynamics allows the formation of a metamor- 

714 phic complex composed of juxtaposed of rocks with different metamorphic 

715 evolutions, likely with a homogeneous lithostratigraphic setting, accreted in 

716 the mantle wedge and exhumed during active subduction. To evaluate the 
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717 compatibility between natural data and the tliermo-mechanical configuration 



718 (generated by alternative models, such as micro-continent collision (Rosen- 



719 jjaum fc Lister , 2005) and intra-continental subduction (Thoni, 2006 Janak 



et a/., 2009), numerical simulations of these models could be performed in 



721 the future using the same approach. 
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FIGURE CAPTIONS 

Fig. 1: Locations from which natural data were obtained within the Aus- 
if roalpine Domain upon the structural map of the Alps, redrawn after Schmid| 
eFal] ( pOOil ), dThonif [2006| and [CasteUi et d] ( |2007| . Western Alps (a). 



^efs. 1: iCasteUil (|1991j); 2: [Vuichard fc Ballevre| (|1988j); 3: [Lardeaux fc 



Spalla (1991); 4: CasteUi & Rubatto (2002); 5: Lardeaux et al. (1982); 6: 



Gosso et al. (1982); 7: Hy (1984); 8: Oberhansh et al. (1985); 9: Compagnoni 



1303 (1977); 10: Zucah et al. (2002); 11: Delleani et al. (2010); 12: Rubatto et al. 



J p99| ); 13: [Giorgetti et ol] ( [2000| ) ; 14: [Tropper fc Essene] ( p002| ) ; 15: [Tropper 



et al. (1999); 16: Williams & Compagnoni (1983); 17: Spalla et al. (1997) 



1306 and Gosso et al. (2010); 18: Ruffet et al. (1995); 19: Zucah et al. (2004) and 



Zucah k Spalla (2011); 20: Andreoh et al. (1976); 21: Reinsch (1979); 22: 



Desmons k Ghent (1977); 23: Pognante (1989); 24: Pognante et al. (1987); 



1309 25: Rebay k Messiga (2007) and Lardeaux k Spalla (1991); 26: Pognante 



1310 (1989); 27: Williams k Compagnoni (1983); 28: Pognante et al. (1987); 29: 



Lardeaux et al. (1983); 30: Pognante et al. (1987); 31: Pognante (1989); 



1312 32: Rolfo et al. (2004); 33: Hellwig (2003); 34: Roda k Zucah (2008); 35: 



^ienast k Nicot] p97l| ) ; 36: |Hopfer| ( |1995| ) ; 37: |Kienast| p983| ); 38: |Dal Piaz 



et al. (1983); 39: Pennacchioni (1988); 40: Scambelluri et al. (1998); 41: 



Dal Piaz et al. (2003); 42: Cortiana et al. (1998); 43 and 44: Rubatto et al. 



1316 (2011). Central Alps (b). Refs. 1: Gazzola et al. (2000); 2: SpaUa et al. 



J p95| ); 3: [Tomaschek k Blumel| ( |1998| ); 4: [Handy et al\ ( |1996| ); 5: [Gunth 



k Liniger (1989); 6: Konzett k Hoinkes (1996); 7: Habler et al. (2001); 8: 



Solva et al. (2001); 9: Spalla (1993); 10: Hoinkes et al. (1991); 11: Zanchetta 



1320 (2007); 12: Habler et al. (2006). Eastern Alps (c). Refs. 1: KoUer (2003); 2: 
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1321 ILinnerl ((19991); 3 and 4 IFaryad & Hoinkesl ([20031); 5: IGaidies aZ.I (|2008aP; 



Gaidies et al. 


(2008b 


); 


Stiiwe & Powell 


(1995) 



Gregurek et al. (1997); 15: Tenczer et al. (2006); 16: Janak et al. (2009); 17: 



Putis et al. (2002); 18: Neubauer et al. (1999); 19: Torok (1999) and Torok 



1326 (2003); 20: Tropper et al. (2001). 



1327 Fig. 2: Cross-sections for Eastern (AA'), Central (BB') and Western (CC) 



1328 ^Ips redrawn after Polino et al. (1990), Dal Piaz et al. (2001) and Schus- 



ter et al. (2004). In the Central Alps cross section the lower austroalpine 



1330 comprises slices of basement units of the Err and Bernina systems, interlay- 



1331 ered with thin slivers of Mesozoic sediments (e.g. Froitzheim et al, 1994), 

1332 giving rise to a complex setting characterised by the repetition of cover and 

1333 basement lithostratigraphic units. 



1334 Fig. 3: P/T ratios and peak ages obtained from the collection of natural 

1335 data (Tables 1 and 2). 



1336 Fig. 4: Peak P-T estimates obtained from the natural data collected for the 

1337 Western (a). Central (b) and Eastern (c) Austroalpine; the numbers refer to 

1338 locations marked in Fig. 1. The metamorphic fades fields shown are after 



1339 Ernst & Liou (2008). GS: greenschist facies; BS: blueschist fades; EP-A: 

1340 epidote-amphibolite facies; A: amphibolite facies; A-E: amphibole-eclogite 

1341 facies; EP-E: zoisite-eclogite facies; Lws-E: lawsonite-eclogite facies; GR: 

1342 granulite facies; HGR: kyanite-granulite facies. 



Fig. 5: Four evolutionary stages of the model: (a) 10 Ma, (b) 20 Ma, (c) 
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1344 35 Ma and (d) 60 Ma. In the insets, the streamhnes (blue hnes) obtained 

1345 in the corner area are reported. 1) Upper continental crust, 2) lower conti- 

1346 nental crust, 3) upper oceanic crust, 4) lower oceanic crust, 5) sediments, 6) 

1347 exhumed upper continental crust, 7) exhumed lower continental crust, and 

1348 8) continental lithospheric mantle. 

1349 Fig. 6: An example of a P-T-t path for a lower continental marker at four 

1350 time increments; the red circle indicates the initial position of the marker, 

1351 and the black circle indicates the marker location at a particular subduction 



1352 age. The metamorphic facies fields are from Ernst & Liou (2008). GS: 

1353 greenschist facies; BS: blueschist facies; EP-A: epidote-amphibolite facies; 

1354 A: amphibolite facies; A-E: amphibole-eclogite facies; EP-E: zoisite-eclogite 

1355 facies; Lws-E: lawsonite-eclogite facies; GR: granulite facies; HGR: kyanite- 

1356 granulite facies. 

1357 Fig. 7: Predicted peak P and T conditions recorded by the exhumed mark- 

1358 ers (see the legend for the lithologic affinities) in terms of PmaxT (a) and 

1359 PTmax (b). (c): thermal gradients in the subduction zone for seven refer- 

1360 ence depths obtained from a comparison of the dynamics evolution of the 

1361 isotherm configuration to the thermal gradients obtained by the burial paths 

1362 of two markers (blue = upper continental crust; green = lower continental 

1363 crust), (d): The P/T ratio obtained based on P/T ratio obtained by the 

1364 mean PmaxT conditions recorded by the exhumed markers for five different 

1365 time steps, a): Upper continental crust and b): lower continental crust. 

1366 Fig. 8: Comparison of the natural P-T peak assemblages coming from 
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1367 the Austroalpine of the Western Alps (blue dots) and simulated PmaxT and 

1368 PTyyidx data (grey dots) for slices ablated from upper (a-c) and lower (b-d) 

1369 continental crust, e) P-T peak assemblages obtained for the peak age of 

1370 SLZ (blue dots) and f) for the External Slices (blue dots) compared to the 

1371 natural estimates (shaded area). The orange dots represent the simulated 

1372 P-T data at the end of the subduction period (30 Ma), which is interpreted 

1373 as the exhumation age. 



1374 Fig. 9: Comparison of the simulated and natural P-T-t paths inferred for the 

1375 Austroalpine. See caption of Fig. 1 for the reference numbers. Pre-Alpine 



( 


1991 


); 25: 


Lardeaux 



& Spalla (1991) and Rebay & Spalla (2001); 34: Kienast & Nicot (1971), 



Gardien et al. (1994) and Roda & Zucah (2008); 01: Gazzola et al. (2000) 



1379 and 19.2: Schuster et al. (2001). 



1380 Fig. 10: Comparison of the natural P-T peak assemblages coming from 

1381 the Austroalpine of the Central Alps (red dots) and simulated Pmax^ and 

1382 PTfYiax data (grey dots) for slices ablated from upper (a-c) and lower (b-d) 

1383 continental crust, e): P-T peak assemblages obtained for the peak age of 

1384 the Central Austroalpine (red dots) compared with the natural assemblages 

1385 (shaded area). The orange dots represent the P-T data for the period from 

1386 50 to 35 Ma, which is interpreted as the exhumation age. 



1387 Fig. 11: Comparison of the natural P-T peak assemblages coming from 

1388 the Austroalpine of the Eastern Alps (green dots) and simulated Pmax^ and 

1389 PT^aax data (grey dots) for slices ablated from upper (a-c) and lower (b-d) 



60 



1390 continental crust, e) and f): P-T peak conditions predicted after 30-40 Ma 

1391 from the beginning of subduction obtained for the two different starting ages 

1392 of the simulation: 130 Ma (e) and 100 Ma (f) compared to natural estimates 

1393 (shaded area). The red and orange dots represent the P-T data at two 

1394 times from the start of the subduction and are interpreted as two natural 

1395 exhumation ages (80 Ma red; 50 Ma orange). 
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Table 1: Climax ages for the Austroalpine of the Western Alps: EMC Eclogite Micaschists 
Complex; GM Gneiss Minuti; II DK II Zona Dioritico-Kinzigitica; RCTS Rocca Canavese 
Thrust Sheet; CH Chatillon; DB Dent-Blanchc; VP ValpcUinc; RS Roisan; EL Etirol- 
Levaz; EM Emilius; GL Glacier-Rafray; PK Pillonet Klippc. Keys correspond to the 
references and locations on Fig. la. 
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Valpclline 


VP 


37 


Etirol-Lcvaz 


EL 


39 


Mt. Emilius 


EM 


41 


Glacier-Rafray 


GL 


42 


Ayas 


PK 


43 


Quincinetto 


EMC 


44 


Brosso 


EMC 



65.0±5.0 
71.2±3.2 
118.0±2.3 
129.0±15.0 
114. Oil. 
69-90 



65.0±3.0 
68.0±7.0 
>80 
70-90 
>69.4±0.7 
65.0±3.0 
> 45-60 
> 60-85 
36.0±4.0 
41.0±0.3 
46. Oil. 5 
46. Oil. 5 
45.0i0.7 
47.0i0.9 
40.0i0.5 
49.0i0.5 
45.0i0.4 
74. Oil. 
78.5i0.9 
76.8i0.9 



Rb/Sr 
Ar/Ar 
Rb/Sr 
Rb/Sr 
Rb/Sr 



Ar/Ar 
Rb/Sr 
Ar/Ar 
zircon 
Rb/Sr 
K/Ar 
K/Ar 
Ar/Ar 
Rb/Sr 
Rb/Sr 
Rb/Sr 
Rb/Sr 
Rb/Sr 
Rb/Sr 
Rb/Sr 
Rb/Sr 



Rubatto et al. 



Hy 



Oberhansli et al. 




1985 


Oberhansli et al. 




1985 



Hunziker 1974 



Rubatto et al. 


1999 


Cenki-Tok et al. 


12011 


Rubatto et al. 


1999 


Rubatto et al. 


1999 



Stockhert et i 



Hunziker 1974 



Ruffet et al. 


1995 


Rubatto et al. 


1 1999 


Hunziker et al. 


1992 


Hunziker et al. 


1992 


Ayrton et al. 


1982 1 


Cosca et al. i 


1994 



Ayrton et al. 



Ayrton et al. 



Dal Piaz et al. 


2003 


Dal Piaz et al. 


2003 


Dal Piaz et al. 


2003 


Dal Piaz et al. 


2003 


Dal Piaz et al. 


2003 


Cortiana et al. 


1998 


Rubatto et al. 


2011 


Rubatto et al. 


2011 



1982 



1982 



74 



Table 2: Climax ages for the Austroalpine of the Central-Eastern Alps: UA Upper Aus- 
troalpine; LA Lower Austroalpine. Keys correspond to the references and locations on 
Fig. lb (Central Alps) and Fig. Ic (Eastern Alps). 



Key 


Locality 


Nappe 


Climax 


Method 


Refer en ecs 






Met. C. 


age (Ma) 






Central Alps 



Maloja 
Engadina 



6 Moos 

7 Schnccbcrg 

8 Montcncve 

10 Saltaus 

11 Ulfas 

12 Saltaus 



LA/EiT 
LA/Margna 

UA/Schneeberg 
U A/ Schnccbcrg 
UA 
Tcxcl G. 
UA/Texcl G. 
UA 
Tcxcl G. 
UA/Tcxcl G: 



76-89 
60-80 
60-80 
84.5±1.0 
91.5-96.3 
>76.1±4.8 
>83.2±2.9 
143.0±2.0 
84.0±5.0 
>95.0±5.0 
85.2±4.6 



K/Ar 
Rb/Sr 
K/Ar 
Ar/Ar 
Sm/Nd 
Rb/Sr 
Ar/Ar 
Rb/Sr 
U/Pb 
Sm/Nd 
Sm/Nd 



Handy et al. 



Frey et aL 



Frey et aL 



1974 



1974 



Konzett & Hoinkes i 1996 



Habler et al. i 2001 



Hoinkes et al. i 1991 



Zanchetta i 2007 



Solva et al. • 2005 



Solva et al. 


2001 


Solva et al. 


2001 



Habler et al. i 2006 i 



Eastern Alps 



1 


Rcckiier 


LA/Reckner N. 


2 


Schober 


UA/Schobergruppe 


3 


Mur Valley 


UA/Woelz C. 


7 


Kupplerbrunn 


UA/Saualpc 


8 


Saualpc 


UA/Saualpe 


15 


Koralpe 


UA/Koralm C. 


16 


Pohorjc 


UA/Koralm C. 


17 


Sieggraben 


UA/Sicggraben C. 


18 


Sieggraben 


UA/Sicggraben C. 



49.5-51.8 

87-115 
92.4±14.4 
91.1±1.3 
90.0±3.0 
87.5±4.5 
92.0±0.5 
>103.0±14.0 
>136-108 



Ar/Ar 
Sm/Nd 
Sm/Nd 
Sm/Nd 
Sm/Nd 
Sm/Nd 
U/Pb 
U/Pb 
Ar/Ar 



Dallmeyer et al. 



1997 



Schuster & Frank 



Thoni et al. 



Thoni & Miller 



Miller & Thoni 



Janak et al. 



Putis et al. 



1996 



1997 



2000 



Dallmeyer et al. 11996 
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Table 3: Material properties used for the simulations. Reference s: (a) ([Ranalli fc Murpliy[ 
[r987|), (b) dKkbyl [19831) (c) dChopra fc Petersonj [l98T| ) (d) (|Karato fc Wu[|1987D, (e) 
([Dubois fc Dianient[|1997||Best fc Christiansen[|2001[ ), (f) (Haenel et aL [ |1988[ ) ' (g) ( [Gerya 
fc Stockliert[|2002|). 



Materials 


Rheology 






n 


po(feg/m^) 


k(W/mK) 




E(KJ/mol) 


Rcfs. 


Continental crust 


dry granite 


3.47 


10^1 


3.20 


2640 


3.01 


2.50 


123 


a,G,f 


Oceanic crust: upper 










2961 


2.10 


0.40 


260 


b,e,f 


Oceanic crust: lower 


diabase 


1.61 


jq22 


3.40 


2961 


2.10 


0.40 


260 


b,e,f 


Sediments 


IQl^ 








2640 


3.01 


2.50 


123 


e,f 


Dry mantle 


dry dunitc 


5.01 


10=0 


3.41 


3300 


4.15 


0.002 


444 


c,e,f 


Serpentinized mantle 










3000 


4.15 


0.002 


444 




Water 










1000 


0.60 








Atmosphere 


J023 








1.18 


0.03 
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